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ABSTRACT 

Redshifts of a supernova (SN) and gamma-ray burst (GRB) samples are compared 
with the pixel temperatures of the Wilkinson Microwave Anisotropy Probe (WMAP) 
seven- years data, the pixels locations corresponding to the SN and GRB sky coordi- 
nates. We have found a statistically significant correlation of the SN redshifts with 
the WMAP data, the average temperature deviation being +29.9 ± 4.4 /iK for the 
redshifts z ranging from 0.5 to 1.0 and +8.6 ± 1.3 nK for z e (0.0,0.4). The latter 
value accords with the theoretical estimates for the distortion of the cosmic microwave 
background due to the integrated Sachs- Wolfe effect, whereas the larger anomaly for 
higher redshifts should be studied in more detail in the future. 
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1 INTRODUCTION 



During the past few years there has been a notice- 
able increase in publications discussing different as- 



crowave background (CMB) radiation (|Ferreira et al. 


1998; 


Wu et al.ll200]J; Verde et alj2001 : Komatsu et al.ll2003 


: Park 


2004 


: Vielva et al. 12004 De Troia et al. 20071: Vanchurin 


2009 


: Bernui & Reboucasl I201QL iPogosvan et all l201ll: 


Ashoorioona et all 201lf). The non-Gaussian behaviour seen 



other known anomalies are relat ed to the so-called "hot" 
and "cold" spots of t he CMB (lLarson fc Wandeltl |2004 
l2005l : IVielva et alll2007l : ICopi et al.ll2010l : lAvaita et alj201Ch 
which are regions of increased or reduced temperatures com- 
pared to the average temperature of the CMB map. For 
example, the temperature of one of the most prominent 
and broadly discuss ed non-Gaussian features known as the 
extreme cold spot (Vielva et al 20041: iMukheriee fc Wangl 



on different angular scales and with a variety of amplitudes 
imply the existence of some large-scale anisotropics in the 
distribution of matter in the Universe. 

One of such anisotropics is a so far unexplained 
alignment of low-order multipoles (I — 2 to 5) in 
the representation of CMB as a power series expan- 
sion with spherical harmonics. This anomaly is sometimes 
referr e d to as "the cosmic axis o f evil" (|Spergel et all 



20031: |Pe Oliveira-Costa et all 12004 ISchwarz et all 12004 



Land fc Magueiiol 120051: iMagueiio fc Sorkinl 120071 : see also 
the review by iRakic fc Schwarj 120071 ). Furthermore, re- 
cently it was found that the axes of rotation of most 
galaxies from the Sloan Digital Sky Survey (SPS S) ap- 
pear to line up with the axis of evil (|Longdl201lh . The 
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|2004 iMartmez et al.ll2006l : ICruz et"ai1l2007l : IVielvall2010l ) 

approximately 70 fiK lower than the average CMB temper- 
ature. The possibilities of such a large anomaly appearing 
within the framework of the sta ndard A-cold -dark-matter 
(ACDM) model are discussed by llnouej (|2012h . who argues 
that there might exist a large-scale void along the line of 
view towards the extreme cold spot. The existence of a 
void of 100 Mpc radius a long this line a t 0.5 < z < 0.9 
was recently ruled out by iGranett et all l|2010t ). who used 
galaxy counts in photometric redshift bins. However, the 
void might be at larger distances. The non-uniform distri- 
bution of the largest CMB spots over the celestial sphere 
results in the north-south ecliptic asym metry in the CMB 
maps (lEriksen et alj|2004 | Bernuill2008l ; see also the discus- 
sion m IRakic fc Schwarzll2007l ) 



A variety of techniques are used for extracting the in- 
formation about CMB anisotropies, such as the determi- 
nation of both angular and planar modulati on parameters 
by u sing the two-point correlation function l|Szapudi et all 
120011 ). the characterisation of deviations from a given (say, 
Gaussian) model using the x 2 -criterion, or the comparison of 
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random Ga ussian maps with observed maps |Liguori et al.l 
I2003l , l2007l ). 

Also of importance is the interpretation of the devia- 
tions from a Gaussian distribution found in the CMB power 
spectrum. One of the possible explanations of these devia- 
tions is based on the assumption that they might be due to 
contamination of CMB by foreground sources. The induced 
deviations caused by the strongest of these sources are re- 
moved by the foreground reduction procedure, although the 
completeness of this procedure is still under scrutin y (see 
iNaselskv et~alll2003l ; lThenll2006l ; iBennett et ai1l201ll ). 

It is wo rt hwhil e mentioning th e Sachs-Wolfe 
llSachs fc Wolfd Il967l ). Ree s-Sciama dRees fc Sciamal 
1 19681 ) and Sunvaev-Zeldovich (|Sunvaev fc Zeldovichl Il970l ) 
effects. The first two alter the energy of CMB photons 
when they traverse gravitational fields of different strengths 
corresponding to either voids or conglomerations of matter. 
In ACDM models, the "late-time " Sachs-Wolfe effect works 
mostly on large angular scales (|Bartolo et alj [2010). The 
third is related to the scattering of CMB photons from high 
energy electrons in hot coronae of galaxy clusters. 

The distinction between the proper and induced CMB 
anomalies can be made by comparing the WMAP data with 
patterns created by other sources on the sky. For example, 
IVerkhodanov et al.l (|2010l ) have found statistically signifi- 
cant correlations between CMB map temperatures and the 
distribution of gamma-ray bursts (GRB). Significant corre- 
lations of WMAP data with galaxy s amples on large angu- 
lar sca les of 4° to 10° were found bv lFosalba fc Gastanagal 
J2004I) for the A u tomat e d Plate Mea s uring galaxy survey, 
by iFosalba et all (l2003h: ICabre et all ||2006| ) for the Sloan 
Digital Sky Survey. bv lGoto et al.l (|2012l ) for the wide-field 
infrared survey of galaxies and by other authors. 

Since there exists a broad variety of opinions with re- 
spect to the nature of the anomalies seen in the WMAP 
data, further statistical studies of these data by using alter- 
native techniques are desirable. In particular, it is important 
to continue with cross-correlations of CMB maps with spa- 
tial distributions of different high-redshift objects, such as 
galaxies, quasars or GRBs. 

Most of the cross-correlations are based on two- 
dimensional distributions of extragalactic sources. Here we 
shall explore the third coordinate (redshift z) for checking 
whether some of the CMB fluctuations might be caused by 
distant concentrations of matter. We shall use supernovae 
(SNe) and GRB samples as the test objects that trace mat- 
ter in the universe. 



2 CMB TEMPERATURE AS A FUNCTION OF 
SN REDSHIFTS 

For our analysis we ha v e ta ken the list of su- 
pernovae (|Bartunov et al.l 120071 ) compiled by the 
Sternberng Astronomical Institute (Moscow) at 
|http: //www, sai .msu.su / sn/ sncat which contains a 
large sample of SNe with known redshifts (5688 SNe for 
September 2011). We have also used sources from the 
GRB list maintained by J. Greiner ICreineij l|2008l ) at Max- 
Planck-Institut fur extraterrestrische Physik (Garching) 
at http://www.mpe.mpg.de/~jcg/grbgen.html. With the 
updates from ICucciara et al.l (|2011l ) and the Swift data 
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Figure 1. Histogram of the WMAP pixel temperatures TgN near 
the supernovae locations as a function of SN redshifts ~sn (for the 
WMAP ILS values Tg^); the points with no error-bars correspond 
to the redshift bins containing a single source. The dashed line 
indicates the slope of the linear regression of the SN temperatures: 
T SN = (+23.7 ±6.0) z SN fj,K 



archive at http://www.swift.ac.uk this list contained 992 
GRBs (as for October 2011). 

The WMAP data products l|The WMAP Group! 
2010) were obtained from the NASA archive at 
[httpT77lambda.gsfc.nasa.gov, and the software for 
processing these data was adopted from the comprehensive 
HEALPix package developed by Jet Propulsion Labo- 
ratory and availa ble at http://healpix.jpl.nasa.gov 
|Gorski et al.l l2005). Since most of the correlations between 
the CMB and extragalactic sources are found on large 
angular scales, we have used here the foreground reduced 
WMAP sky maps with the resolution corresponding to the 
HEALPix parameter A sidc = 128. 

Figure [JJ shows the histogram of the WMAP pixel tem- 
peratures for the SNe locations as a function of SNe redshifts 
for the WMAP Internal Linear Combination (ILC) map with 
the redshift ranging from zero to 1.4. 

In order to avoid contamination from point sources and 
from the Galactic pla ne, we have used the WMA P temper- 
ature analysis mask (|The WMAP Group! |2010| ) which in- 
cludes point sources from external catalogues, as well all as 
the proper WMAP point source catalogue . It also takes int o 
account the galactic foreground emission |Gold et al.ll201l1 ). 
We have taken the temperature analysis mask with higher 
resolution corresponding to A^idc = 512 (Figure [2} and ap- 
plied it to the same higher resolution ILC-map, after which 
we degraded the masked ILC-map to the lower resolution of 
A sidc = 128. 

Additionally, we have restricted our SN and GRB sam- 
ples to the galactic latitudes \b\ > 40° and rejected the 
WMAP pixels with temperature values |T| > 0.5 mK, as 
such pixels are likely to correspond to artifacts or some extra 
point sources. The distribution of SNe over the sky within 
the selected latitude limits is shown in Figure [3] 

As for the GRB data, only part of them contains pho- 
tometric redshifts, so we included in our a nalysis the GRB 
pseudo-redshifts from the catalogue of IPelangeon et al.l 
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Figure 2. WMAP temperature analysis mask applied to the 
WMAP ILC-map prior the computation of the SN and GRB red- 
shift histograms. 




Figure 3. Sky distribution of SNe for the selected range of lati- 
tudes |fe| > 40°. The strip showing high SN densities corresponds 
to the SDSS Supernova Survey covering the equatorial strip of 
approximately 2.5°-width 



(2008) available at |http : / /cosmos . ast . obs-mip . f r| (the 
empirical methods for estimating pseud o-redshifts are de- 
scribed bv lAtteia|[2005l ; lBand et alj|2004h . 

By grouping the WMAP pixel temperatures into red- 
shift bins with the bin size of 0.2 for GRBs and 0.01 for 
SNe, we have built two histograms for examining the pos- 
sible relationship between the SN and GRB redshifts and 
CMB temperatures. 

The histogram in Figure [T] shows the bin-averaged val- 
ues of the WMAP temperatures Tsn for supernovae; the 
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Figure 4. Histogram of the WMAP pixel temperatures Tgrb 
near the GRB locations on the sky as a function of redshifts zq^b 
(for the WMAP ILC values Tgrb); the points with no error- bars 
correspond to the redshift bins containing a single source. 



values Tgrb for gamma-ray bursts are shown in Figure [4] 
The error-bars in these plots indicate the standard errors 
(SE) of the averages Tsn (or Tgrb) in each redshift bin cal- 
culated as SE = v/XJ-LoC^In — Tsn) 2 /n(n — 1), where n is 
the number of SNe in the bin. The points corresponding to 
a single source in the bin are shown having no error-bars. 
The average WMAP pixel temperature for the selected lat- 
itude range, T\ b \ >i0 o — +2.7 ±0.2 fiK, was subtracted from 
the bin-averages (the positive shift of T\ b \ >40 o is due to the 
latitude selection effect, as the zero-average corresponds to 
the whole sky only). Although the data are noisy, one can 
notice a statistically significant deviation of Tsn from the 
average temperature Ti b i >4a o, which is larger for the red- 
shifts ranging from 0.5 to 1.0, where the average of Tsn 
is Tsn 5 ' 1 ° — +29.9 ± 4.4 /iK, and which is smaller for 
lower redshifts z £ (0.0,0.4): Tsn ' ' 4 ' = +8.6 ± 1.3 fiK, 
the average number of SNe in each redshift bin being 52 for 
2 £ (0.0, 0.4) and 6.5 for z £ (0.5, 1.0). For the whole sample 
of 2413 SNe within the selected latitude range, the average 
temperature is also significantly positively biased: Tsn = 
+ 12.0 + 1.6 fjK. Such an excess of Tsn is larger than the cur- 
rent theoretical and observational estimates for the potential 



tamination J Knox et al. 


1998l:ISantos et al.ll2003l:IZahn et all 


120051: ISerra et al. 20081: 


Taburet et alj|2009h. 



The GRB histogram (Figure [4} is quite random and 
shows no temperature excess. For this histogram Tgrb = 
+ 1.3+5.7 fjK, which is statistically insignificant and which is 
what one would expect [the average number of GRBs in each 
redshift bin was 6.4 for the redshift range z £ (0,4)]. One 
could attribute the difference between the two histograms 
to the differences between the GRB and SN distributions 
across the sky: the GRB sample is uniform, whereas the 
SN sample contains fields with high source densities, cor- 
responding to specialised surveys of extragalactic sources. 
Figure [3] shows the sky distribution of the SNe used here 
for calculating the temperature-redshift histogram. There is 
a distinctive strip seen in this map that contains a larger 
number of SNe. It corresponds to the SDSS Supernova Sur- 
vey conducted for a narrow equatorial zone of ~ 2.5°-width 
l|Frieman e t al. 2008). In order to check whether the excess 
of the WMAP pixel temperatures is due to this SN subsam- 
ple, we have calculated a separate TsN/zsN-histogram con- 
taining only the SNe with \S\ < 1.5° (shown in Figure [5}, 
and we have found that this histogram reveals a statistically- 
significant positive shift of the average temperature for low 
redshifts Tg^ 0,0 ' 4 = +11.1 + 1.4 fj,K, whereas for higher red- 
shifts the shift is not significant: Tsn 8 ' 1 ° = +3.6 ± 4.7 pK, 
which indicates that the found temperature shift for higher 
redshifts comes from the rest of the SN sample, which is 
distributed more uniformly across the sky than the SDSS 
Supernova Syrvey. 

In order to distinguish between possible sources of the 
detected anomaly we have calculated similar histograms for 
three different WMAP frequency bands, separately for the 
Q, V and W bands (keeping in mind that, for instance, the 
Sunyaev-Zeldovich effect is energy-dependent). These three 
histograms are similar to each other and to the histogram 
shown in Figure ffj their linear regression slopes varying only 
slightly: +27.4 ± 6.9, +35.6 ± 7.6 and +32.2 ± 7.8 for the Q, 
V and W bands, respectively. This means that a larger sam- 
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3 DISCUSSION 
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Figure 5. Histogram TgN/ z SN f° r the supernova subset within 
the equatorial latitudes ±1.5° seen as a high-SN-density strip in 
Figure [3] 
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Figure 6. Histogram Tsn/zsn for the randomised sample of SNe 
(to be compared with the histogram shown in Figure [3] 



pie of SNe (which might become available in the future) 
is needed to get a more clear picture about the nature of 
the detected correlation. A similar band-independent corre- 
lation between WMAP temperatures and large-s c ale st ruc- 
tures traced by galaxies was found bv lGoto et all |2012l ) for 
the Wide-field Infrared Survey Explorer data, which was at- 
tributed by the authors to the integrated Sachs- Wolfe effect. 

Yet another check ensuring that the seen correlation is 
not due to a systematic error consists in randomising the SN 
positions on the sky. For this purpose, we have used the same 
SN sample with assigned random coordinates, while main- 
taining their redshift-distribution, as well as their higher 
density in the equatorial zone corresponding to the SDSS 
subsample. The histogram for this randomised SN loca- 
tions (Figure [6} does not reveal any significant temperature 

-1.7±1.8 fiK) and high 
+9.3±6.6 fiK) redshift ranges, with the average 
for the whole SN sample being TsN rand = — 0.3±1.6 /iK. So, 
the found anomaly is likely to be real. 



anomalies, for both low (Tsn°' a 

,=(0.5,1.0) 
V 1 SN rand 



In our view, this anomaly can be interpreted in a number 
of different ways, of which we would like to mention the 
following: 

- it might correspond to the secondary CMB anisotropics 
caused by dust associated with clusters and superclusters of 
galaxies at moderate redshifts z ~ 1; 

- it could be related to secondary and tertiary ionisation 
of the foreground media during the post-reionisation epoch; 

- it could be caused by large-scale inhomogeneities in 
matter distribution, such as voids, filaments or walls, which 
could alter the energies of the primordial CMB photons 
through the integrated Sachs- Wolfe and Rees-Sciama effects. 

These effects are broadly discussed in the literature, so 
let us just briefly outline how they contribute to the CMB 
anisotropies. 

Dust contribution. The SN and GRB host galaxies con- 
tain dust whose thermal radiation contributes to the CMB. 
It is known th at dust in sp i ral galaxies has a wide range o f 
temperatures (|RowarJll99>j iFranceschini fc Andrean3ll995t ) 
which makes it difficult to separate CMB anisotropies 
caused by dust from the primordial CMB anisotropies. Re- 
cent Hershel submillimiter observations have shown the 
presence of large dust masses in high-redshift galaxies 
(ISantini et al. I l2010h . Similarly, large dust masses have been 
measured in loca l universe galaxies (I Galliano et all | 2005| ; 
Bendo et al.ll2006l: Icialametz et alj|2009l : iGrossi et al J 1201(3 : 



O'Halloran et al.ll2010h . So. indeed, there is a possibility that 



dust from galaxies at all redshifts contributes to the CMB 
temperature fluctuations, which, in turn, would result in a 
positive correlation between C MB tempe r ature s and over- 
densities of galaxies - see, e.g., iHo et all ((2008). The pos- 
sible impact of dusty galaxy clusters o n cosmological pa - 
rameter estimates wa s als o discussed b y Serra et all (|2008h . 
iDunklev et"afl (|201ll ) and lMillea et all (|2012r i. 

Ionised gas. Another possible source of secondary 
anisotropies in the CMB power spectrum is the Thomson 
scattering on free electrons in regions of ionised gas with 
bulk pe culiar velocities produced at the epoc hs of reion- 
isation |McQuinn et al. | 120051 : iMcQuinnl l2010l l and post- 
reionisation. Local variations in the ionised gas increase 
temperature anisotropies through the kinetic Sunyaev- 
Zeldovich effect, the primordial temperature perturbations 
being damped as e~ T , where r is the Thomson scattering 
optical depth. 

The WMAP best-fitting model for the large-scale po- 
larisation an isotropy favours a mean redshift of reionis ation 
of 10.4±1. 4 jBennett et ap2003l : lKomatsu et al.ll2009l ') with 
r w 0.08 l|Dunklev et al-lhoOSh . The observed spectra of 
quasars indicate a rapid change in the ionising background 
at the redshift z ~ 6, so that reionisation is completed at 
that epoch. By contrast, there are many uncertainties re- 
garding the redshift corresponding to the beginning of reion- 
isation. The polarisation anisotropies models are based on 
a number of assumptions regarding the density distribution 
of the baryons, quasar physics, etc., much of which are un- 
known. If reionisation occurred gradually ove r a period of 
time, then it started as long ago as z = 30 (|Kogut et al.l 
l2003h. so that the values of r = 0.08 or even r = 0.17 
( Kogut et al.ll2003l : ISpergel et all 120031 ) is likely to be un- 
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derestimated. As was shown bv lGnedir] (|2004t ). one can con- 
struct a posteriori models that produce a large value of the 
Thomson scattering optical depth (up to r = 0.2) consistent 
with the WMAP measurements. Moreover, the reionisation 
history could be much more complex than is usually thought 
dBarkana fc Loebll200ll ; lMiralda-Escude1l2003l; lBrommll2004l ; 



McQuinn et al.ll2007t). Accord i ng to Wvithe fc Loebl l|2003f ); 
Cenl (|2003l l; lHiman fc Holderl 120031 ); iBromml (|2004h . there 



might have been more than one reionisation period, with 
the first reionisation occurring at z ~ 20 due to massive 
population III stars, followed by a period of partial recom- 
bination, and, finally, the reionisation at z ~ 6. Giant blobs 
of gas are known to surround galaxies and galaxy clusters 
llFrancis et al.lll996l ; ISteidel et aifeoool ; iMatsuda et"al]|2004l ; 
lHaves et al.ll201lT l. Some of partially ionised blobs remaining 
after the end of the reionisation epoch at z = 6 could screen 
and dampen the CMB fluctuations along certain directions 
at much lower redshifts (this might be yet another possible 
explanation of the extreme Cold Spot anomaly). 

Besides the afore-mentioned damping of the primordial 
fluctuations, there exists an effect caused by inhomogeneous 
reionisation which leads to the generation of new temper- 
ature fluctuations on small an gular scales (|Babich fc Loebl 
l200fj , 120071 ; iHolder et all 120071 ') due to the bulk motion o f 
electrons in overdense regions (|Ostriker fc Vishniad fl98^) . 
The patchiness of reionisation remains theoretically uncer- 
tain and observationally unconstrained, whereas its contri- 
bution to the CMB fluctuations is likely to be comparable 
to, or larger than, that of the Thomson scattering effect. 

The upper limits of ~ 6 /iK on the CMB anisotropics 
due to the secondary thermal and kinetic Suny aev- 
Zeldovich effects were estimated bv lLeuker et al l (|20ld ) and 
IShirokoff et al.l (|201ll ). These are smaller than the amplitude 
of the anomaly we have reported here. 

Large-scale structure. Both CMB anisotropy and large- 
scale structures of the universe (mapped by galaxy surveys) 
are related to the genesis and evolution of primordial pertur- 
bations and to the cosmological model and its parameters. 
However, due to the existence of large-scale structures in the 
universe (high-density knots, filaments, walls etc.) there are 
gravitational potential wells all over the sky which distort 
the CMB power spectrum on large angular scales through 
the integrated Sachs-Wolfe and Rees-Sciama effects. The 
existence of correlations betwee n WMAP maps and large 
scale structures is well-known dCrittenden fc Turok 19961; 
Kinkhabwala fc Kaminkowskil Il999l; Boughn fc Crittendedl 



20041 ; ICabre et al.l 120061 : IStefanescu fc Popall2010h , so the 
correlation whi ch we are reporting here is likely to be of 
a similar kind. iGranett et al.l (|2008T ) measured the ampli- 
tude of the temperature deviations due to the integrated 
Sachs- Wolfe effect for galaxy superclusters in SDSS, which 
was found to be of of the order of +9.6 ± 2.2 fiK. This is 
in a good agreement with the reported here anomaly for 
low redshifts: TgN°'° 4 ' = +8.6 ± 1.3 fj,K, whereas the larger 
anomaly for higher redshifts, T^ 5 ' 1 ' 0) = +29.9 ± 4.4 /j,K, 
exceeds the theoretical amplitudes for the integrated Sachs- 
Wolfe effect. This should be studied in more detail in the 
future. 

As we have seen, many factors of the foreground con- 
tamination could distort the primordial angular power spec- 
trum of CMB fluctuations. All these effects create a rather 



thick patchy screen between the epoch of last-scattering and 
the current epoch. This screen might be responsible for a 
great deal of the observed CMB anisotropy. 



4 CONCLUSIONS 

The correlation between the SNe redshifts and WMAP tem- 
peratures indicates that the WMAP data are still contami- 
nated by extragalactic foreground, despite being thoroughly 
cleaned from known foreground sources. Thus, further stud- 
ies are important in order to clarify the question concerning 
the fraction of the primordial fluctuations remaining in the 
CMB maps, as the knowledge of the power spectrum of these 
fluctuations is pivotal for the calculation of the parameters 
of the standard cosmological model. 
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